Normal modes of a quasi-one-dimensional multi-chain complex plasma 



G. Piacente 1 , F. M. Peeters 1 and J. J. Betouras 1 ' 2 
1 Department of Physics, University of Antwerp (Campus Drie Eiken), B-2610 Antwerpen, Belgium 
2 Max Planck Institute for the Physics of Complex Systems, 
Noethnitzer Strasse 38, Dresden 01187, Germany 
(Dated: February 2, 2008) 

We studied equally charged particles, suspended in a complex plasma, which move in a plane 
and interact with a screened Coulomb potential (Yukawa type) and with an additional external 
confining parabolic potential in one direction, that makes the system quasi-one-dimensional (QID). 
The normal modes of the system are studied in the presence of dissipation. We also investigated 
how a perpendicular magnetic field couples the phonon modes with each other. Two different ways 
of exciting the normal modes are discussed: f ) a uniform excitation of the QID lattice, and 2) a 
local forced excitation of the system in which one particle is driven by e.g. a laser. Our results are in 
very good agreement with recent experimental findings on a finite single chain system (Phys. Rev. 
Lett. 91, 255003 (2003)). Predictions are made for the normal modes of multi-chain structures in 
the presence of damping. 



I. INTRODUCTION 

Since the first experimental observations of the for- 
mation of Coulomb quasi-lattices [ll involving highly 
charged dust particulate in 1994 3, the research 

field of complex plasmas has seen a sustained growth. 
Complex plasmas consist of micrometer-sized ("dust") 
particles immersed in a gaseous plasma background. 
Dust particles typically acquire a negative charge of sev- 
eral thousand elementary charges, so they interact with 
each other through their strong electrostatic repulsion. 
In the experiment, particles are trapped in a horizontal 
layer by a shallow parabolic well, due to two electrodes, 
and can be suspended in a sheath above the electrodes, 
where the gravity force is balanced by the electrostatic 
force. When the electrostatic energy of neighboring par- 
ticles exceeds the thermal energy by an amount T, the 
particles arrange themselves in regular, solid-like struc- 
tures, i.e. Wigner crystals Q- The mutual Coulomb re- 
pulsion of the dust grains is partly screened by the polar- 
ization of the surrounding plasma particles, mostly by the 
gas ions which represent the major species in the sheath. 
Therefore, the average interparticle potential can be well 
represented by a Yukawa (Debye-H-itckel) potential 6] . 

Complex plasmas provide a new system for the study 
of classical crystalline and liquid dynamics and the melt- 
ing processes. For particle size of the order of /mi, the 
dynamical behavior can be monitored directly with the 
use of optical microscopes pj. 

In the present work we study thoroughly the normal 
modes of a classical quasi-one-dimensional (QID) multi- 
chain complex plasma. Such a QID system was recently 
realized experimentally by giving a proper shape to the 
electrodes 0, @- Experimentally many other quasi- 
one-dimensional or strictly onc-dimcnsional systems have 
been realized over the years. Colloidal particles, sus- 
pended in aqueous solution, can be trapped in a poten- 
tial well created by two counterpropagating laser beams 
which form a one dimensional coupled array Q. A 
Coulomb chain confined in a storage ring [Tol | , as well as 



ordered electrons on microchannels filled by liquid He- 
lium are other examples of QID classical systems. 
The latter is one of the candidates to be used for quan- 
tum computing 

El 

If3 | . On the atomic scale a chain- like 
system can be found in compounds such as Hg 3 _ s AsFe 
h4| and in low dimensional systems formed on surfaces 
[15J . A one dimensional chain of gas atoms adsorbed by 
carbon nanotubes has been produced in laboratory and 
its phonon spectrum has been calculated theoretically as- 
suming a Lennard- Jones interaction potential [la. fl7| . 

The classical model, we propose in the present paper, 
reveals a non trivial phase diagram at zero temperature 
and allows to calculate dispersion relations for the normal 
modes, which can be directly investigated experimentally. 
Several generic aspects of the model were investigated re- 
cently [13. Here, the main focus is on the normal modes 
of the system and how they depend on different physical 
situations, e.g. frictional forces and the way they are ex- 
cited. We will make connection with recent experimental 
works 0,[ig|. 

The structure of the paper is as follows: we first give 
in Sec. II an overview of the model, stressing the ground 
state properties. In Sec. Ill we present what has been 
known so far for the normal modes of the system, adding 
more results for clarifications, and then we turn to the 
new results with respect to the presence of the gas drag 
(friction) with or without an applied magnetic field. Fi- 
nally, we turn our attention to the forced oscillations in 
Sec. V, where we discuss recent experiments on normal 
modes in single chain systems. Prediction for multi-chain 
configurations are presented in Sec VI. Finally we con- 
clude in Sec. VII. 



II. MODEL AND GENERAL PROPERTIES 

We consider a system of equally charged particles with 
coordinates fi = (xi,yi) moving in a plane and inter- 
acting with each other through a Yukawa-type potential 
(the screening length A is an external parameter which is 
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FIG. 1: The T = structural phase diagram as a function of 
the inverse screening length k and the density n e . The plotted 
quantities are dimensionless, as it will be for all the figures in 
the paper. 



measured in the experiment H3) an d are confined by a 
parabolic potential in the y-direction. The dimensionless 
Hamiltonian of the system is given by: 



H' 



E 



exp{— — r'j\) 
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where H' = H/Eq, k = ro/X and r 1 — ~r*/ro, 
with ro = (2g 2 /roewg) 1 / 3 as unit of length and E Q — 
(mwog 4 /2e 2 ) 1 / 3 as unit of energy; m and g are the mass 
and the charge of the particles, respectively, e is the di- 
electric constant of the medium the particles are moving 
in and coq measures the strength of the confining poten- 
tial. The dimensionless time is defined as t' = u>ot. Fi- 
nally, it is possible to define a dimensionless temperature 
as T> = T/T with T = E /k B = {muj 2 q 4 /2e 2 ) 1 ^ 3 k B 1 . 

In our previous work |18| we investigated the ground 
state and the melting of this Q1D system. We summa- 
rize here the main results, which we will need in the next 
sections. At T = the particles crystallize in a chain- 
like crystal structure, with a linear density equally dis- 
tributed among the chains. In the case of multiple chains, 
if a is the separation between two neighboring particles 
in the same chain, the chains are staggered by a/2 in 
the x-direction, because this arrangement minimizes the 
electrostatic repulsion. The results for the ground state 
configuration are summarized in the phase diagram de- 
picted in Fig. 1. 

The dimensionless linear density is defined as n e = 
uro/a, where v is the number of chains. For low den- 
sities the particles crystallize in a single chain; with in- 
creasing density a continuous transition ("zig-zag" p^l 
occurs and the single chain splits into two chains. Fur- 
ther increasing the density we found the remarkable be- 



havior that the four-chain structure is stabilized before 
the three-chain structure. The 2 — ► 4 chain transition oc- 
curs through a "zig-zag" transition of each of the chains 
accompanied by a shift of a/4 along x. This 4 chain 
configuration has a relatively small stability range after 
which it transits to a 3 chain configuration through a dis- 
continuous, i.e. first order phase transition. For higher 
values of the density, the 4 chain configuration attains 
again the lowest energy. A further increase of n e will 
lead to more chains, that is 6, 7 and so on. The struc- 
tural transitions are discontinuous, i.e. first order, except 
for the 1 — > 2 transition. 

Raising the temperature the ordered structure melts. 
We studied in detail the melting for such a kind of system 
in Ref. [18]. Due to the anisotropy in the two directions, 
a different behavior of the system in x and y directions 
was found. Two different melting temperatures, T x and 
T y , can be assigned. The main features of the transi- 
tion from the ordered state to the liquid state are: i) a 
reentrant behavior as a function of density; ii) a region 
in density for which the system melts first in the uncon- 
fined direction and then in the confined one: this regime 
resembles the findings of Ref. [5l| in the floating solid 
regime; iii) the reentrant melting occurs near the struc- 
tural transition points. For the nearly Coulomb limit 
(k=0.01) there is no evidence of anisotropic melting, that 
is the system behaves more isotropically. Furthermore, 
the Coulomb system has a melting temperature which is 
on average 15-20 % higher than for the screened Coulomb 
inter-particle interaction with k = 1. On the other hand 
for higher value of K, the system behaves more anisotrop- 
ically and the difference between T x and T y is enhanced. 
In the case k = 3 the melting temperature is on average 
10-15 % lower than for the case k = 1. 



III. NORMAL MODES 

In the present paper we are interested in the normal 
modes of the chain structures, and in particular how 
these modes are modified in the presence of frictional 
forces. We will consider: 1) the uniformly damped mo- 
tion of a normal mode, and 2) the damped propagation 
of a local forced oscillation of a single particle. For these 
purposes we review shortly the normal modes in the ab- 
sence of friction. 



A. Dispersion relations in the absence of friction 

In the absence of drag due to the ion gas and exploiting 
the standard harmonic approximation, the equations of 
motion for small oscillations about the lattice equilibrium 
positions in dimensionless units are in the single chain 
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where U = exp(— — ^ |)/|r^ — r'A is the interparticle 
interaction potential. Considering the translational in- 
variancc of the system along the x direction, we search 
for solutions in the form 



{ x 'm V'n) ex P [Okna - wt)\, 



which results into 
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where D a pij is the dynamical matrix, that is the ma- 
trix of the second derivatives of the Yukawa potential, 
calculated at the equilibrium configuration; Qpj is the 
displacement of particle j from its equilibrium position 
in the j3 direction; (a,0) = (x,y), S a /3,ij, % and Spy are 
unit matrices, in particular 5p y takes into account the 
effect of the confining potential. All the frequencies are 
measured in unit of loq . 

The number of chains determine the number of parti- 
cles in each unit cell and therefore the number of degrees 
of freedom per unit cell. So if I is the number of chains 
there will be 21 branches for the normal mode disper- 
sion curves. Note that for ordinary bidimensional crys- 
tals there are 2 acoustical branches and 2p — 2 optical 
branches p3| , if p is the number of atomic species in the 
unit cell. 

Solving explicitly Eq. (4) for the single chain config- 
uration we obtain that the acoustical and optical eigen- 
frequencies are given respectively by: 



where k is the wave number. Numerical results for the 
dispersion relations were presented in Ref. [18]. 

It is interesting to notice that for the acoustical branch 
the dispersion is positive, that is phase and group veloc- 
ity have the same sign, while for the optical branch the 
(2a) dispersion is negative, i.e. the group velocity is nega- 
tive. Physically the negative dispersion for the single 
chain optical branch can be understood considering that 
the electrostatic repulsion acts oppositely to the force of 
the confining potential and this reduces the oscillation 
frequency with increasing k. 

Another notable feature is the softening of the optical 
branch, accompanied by a hardening of the acoustical 
branch at the values of n e and n where the 1 — > 2 struc- 
tural transition is observed (see Fig. 8 of Ref. [18]), 
which confirms that 1 — > 2 is a continuous transition. 

When a magnetic field is applied in the perpendicular 
direction to the plane the particles are moving in, the 
equations of motions are modified and ylto'c is added to 
the right hand side of Eq. (2a) and —x\iJ c to the right 
hand side of Eq. (2b), where cJ* c = qB/mc is the cy- 
clotron frequency and uj' c = ui c /ujq. It is known p4| that 
in a classical system an external magnetic field does not 
alter the statistical properties of the system and conse- 
quently the structural properties should be insensitive to 
the magnetic field strength. On the other hand the char- 
acter of motion of the particles is altered significantly 
because now the x and y motion are coupled. The spec- 
trum of an infinite bidimensional crystal in a magnetic 
field was obtained in Refs. [25l l2r| . Following Ref. (2(| , 
the dispersion relation for our Q1D system in the pres- 
ence of a perpendicular magnetic field B is obtained from 
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B. Dispersion relations in the presence of friction 

In laboratory experiments of a dusty plasma the parti- 
cles experience a frictional drag due mainly to the back- 
ground neutral gas as well as ions. This drag has a signifi- 
cant effect on the dispersion curves of the normal modes. 
In order to compare experimental data with theory, it 
is necessary to develop a theoretical model in which the 
structure of the crystal as well as damping are included 
as essential elements. This can be easily done by adding 
explicitly the friction term in the equations of motion. 
For the single chain configuration 7'i- should be added 
(5a) to the left hand side of Eq. (2a) and j'y ■ to the left hand 
side of Eq. (2b), where 7' = j/ojr, is the dimensionless 
frictional drag coefficient. Similar equations hold nat- 
urally in the multi-chain structures. The equations of 
motion for the 2 and 3 chain structures are reported for 
completeness in Appendix A. 

Proceeding as before, in this case the eigenfrequencies 

(5b) 

are determined by: 
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For a single chain Eq. (7) gives explicitly: 
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from which we obtain the solutions: 
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The analytical expression for the 2 and 3 chain eigenfre- 
quencies are reported in Appendix B. 

In the limit of small wavenumber k, we find for K/n e 3> 
1 that Eqs (9a) and (9b) reduce respectively to: 
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FIG. 2: Dispersion curves for the normal modes in the pres- 
ence of friction for different values of parameters for a one (a), 
three (c) chain structure. Dependence of the acoustical and 
optical branches in the one chain case on the friction coeffi- 
cient at the critical density where the softening of the optical 
mode is observed is given in (b). 
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The real part of the frequency corresponds to the os- 
cillatory motion while the damping in the time domain 

(10b) ^ s gi ven by the imaginary term £7/2. Note that friction 
also effects value of the frequencies of normal modes. In 
Fig. 2 the normal mode spectra are reported for the dif- 
ferent configurations of the system for different values of 
n e , k and 7'. We used in our calculations values for the 
parameters inferred from the experiment 0. The be- 
havior of the dispersion curves reflects rather closely the 

(11a) case without damping. Some new features should, how- 
ever, be stressed: i) the effect of friction results in general 
into a reduction of the frequencies of vibration both for 
the longitudinal motion and for the transversal one; ii) 
for very small values of the wave number the acousti- 
cal vibrations cannot be excited, which implies that they 

(lib) are overdamped. Such waves can only be excited when 
k > k*(n e , K, 7'); hi) the softening of the optical mode at 
the critical density n* for the transition 1 — > 2 depends 
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FIG. 3: The same as Fig. 2 but now a magnetic field of 
stregth u) c /(j1q = 1 is present. The insets depict the damping 
of magneto-phonon modes. 



on 7', in particular the presence of friction reduces the 
value of n* (see Fig. 2(b)). 

When we include a magnetic field the damped normal 
modes are determined by: 



[(uj 2 - Spy + iju))S a p ti j- 
-D a p.ij + iwuj c £ a pSij]Qpj - 



0. 



(12) 



The corresponding dispersion curve for the single and 
multi-chain structures are reported in Fig. 3. The be- 
havior of the curves resembles the case without damping 



but with an additional shift in frequency due to friction. 
Note that in this case it is no longer possible to obtain 
the phonon frequencies analytically. The anticrossing be- 
tween the two branches in the 1 chain configuration is 
still present (see Fig. 3(b)), as in the case without fric- 
tion (see Fig. 11 of Ref. It is remarkable that the 
cyclotron motion and the friction are coupled and the 
magnetic field introduces a dispersion in the imaginary 
part of uj as well. Now, Im(u>) is no longer constant as a 
function of the wavevector as in the case without a mag- 
netic field. Friction mainly alters the acoustical branches 
of the magneto-phonon modes for small wavevectors. 



IV. FORCED OSCILLATIONS IN A SINGLE 
CHAIN STRUCTURE 

In the experiment of Rcfs. 0, |H |H |H HJ H3 the 

system is set into oscillation by an external driving force 
which acts on the system continuously. The frequency 
of such a forced oscillation is then determined by the 
frequency of the driving force and not by the resonant 
frequencies. This is the effective situation in the exper- 
iments where particle motions are excited by laser ma- 
nipulation, which makes it possible to excite and test 
the dispersion relations of certain types of lattice wave 
HE HE EH ? which are longitudinal waves and, most 
recently also, transverse waves were observed HE! OH. 
Laser light exerts a radiation pressure on the particles 
with a magnitude proportional to the laser intensity 013 . 
In these cases the frequency is purely real since the modes 
are driven. What is observed in the experiments, is that 
as the wave propagates it is spatially damped, which can 
be interpreted in term of a complex wave number [2?) 
k = k r + iki. Following this idea and considering that 
excitations take place when the driving frequency is close 
to the free frequency of the modes, we may neglect to first 
approximation for the single chain structure the external 
force and we looked for particular solutions of the equa- 
tions of motion, in the form: 



(x n ,y' n ) oc exp [i(k r na — cut)] exp (—kina) 



(13) 



as it was done in Ref. for the theoretical calculation 
of the optical branch. This yields for the acoustical and 
optical branch, respectively: 
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FIG. 4: Dispersion relations for a single chain structure as a 
function of the: (a) real part and (b) imaginary part of the 
wavevector. The curves show a strong dependence on N.N., 
the number of nearest neighbors included in the calculations, 
at the edges of the first Brillouin zone. 



This clearly shows that inapplicability of this approach 
to obtain the phonon spectrum in the presence of friction. 

The reason for the divergence of the sums in Eqs. (14a) 
and (14b) is a consequence of the fact that the last ex- 
ponent in Eq. (13) blows up for negative values of n. 
This would suggest that alternatively we should look for 
solutions of the equation of motion in the form: 

{x' n ,y' n ) oc exp [i(k r na — tut)] exp (— ki\n\a), (15) 

i.e. damped waves propagating form the location of ex- 
ternal excitation. But in this case in the imaginary part 
in Eqs. (14a) and (14b) the hyperbolic sine term is re- 
placed by exp (— fci|j|a) and consequently the sum gives 
zero. As a result, Eqs. (14a) and (14b) do not have any 
real solutions for the phonon frequency and this approach 
also fails. 

In order to explain some recent experimental results on 
the transversal modes of a finite one dimensional chain 
, excited by striking one particle with two counterprop- 
agating laser beams such that the effective force acting 
on the particle is Iq sin cut, with Iq the intensity of the 
beam, we have followed another approach. We first con- 
sider a single finite chain of N particles confined in the 
y direction. On one of the particles a time varying force 
is acting. We studied the small displacements from the 
equilibrium configuration of each particle, limiting our- 
selves to first neighbor interactions, which is valid for 
ka > 1. The equations of motion for such a system are: 
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dt'- 
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Requiring the frequency w to be real, both equations 
generate a system of two non linear equations, for k r and 
ki. The results of the calculation are reported in Fig. 
4. This approach has, however, some limits. First of 
all, once the laser acts on a specific particle 0, it is no 
longer possible to consider all the particle as identical, i.e. 
the presence of an external force breaks the symmetry 
of the system which is taken into account in Eq. (13) 
by considering the driven particle as being x n= o = 0. 
Second, the convergence of the series in Eq. (14a) and 
Eq. (14b) is no longer guaranteed. The condition that 
must be fulfilled in order to have a convergent sum is 
kid < K,/n e . As seen from Fig. 4(b), this condition is 
not always satisfied (in the specific case considered in 
the picture the condition for convergence is fc^a < 2). 
Note that the dispersion curves depend strongly on the 
number of terms considered in the sum at the edges of 
the first Brillouin zone (Fig. 4(a)), while it is practically 
independent of the number of neighbors considered in the 
sum in the middle of the first Brillouin zone. The system 
of equations arising from Eqs. (14a) and (14b) is not 
defined when ka — and ka = 7r, because in this case 
the equation for the imaginary part is identically zero. 
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with I — 1, 2, ...,N and F ' y the dimensionless strength 
of the driving force. In order to excite the longitudinal 
vibrations we have considered a force directed along x 
and acting on one of the extremity of the chain, while to 
excite the transversal modes the force acts on the particle 
in the middle of the chain and with y component only, 
as was done experimentally in Ref [7]. Looking for a 
particular solution of Eqs. (15a-b) of the form: 



(x' l ,y' l ) = (Af,Af)e- 



(17) 



we obtained the following set of inhomogeneous linear 
equations for the displacements Ai: 



P,A*i_ x + {uu 2 ac + i~fu ac - 2/3 1 )A?+ 



FnS, 



°',JV 



0. 



2 AI, + «, t 



1 + ijuj opt - 2f3 2 )A v l + 
-foA\ +l -F»5 l „=Q, 



(18a) 



(18b) 



7 



where f3i = n\e~ K / n ' ! (2 + 2n/n e + n 2 /n 2 e ) and /3 2 = 
^3 e -K/n e (i_|_ K /^ e ). The solution to these equations may 
easily be obtained from Kramer's rule |34j : 
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where is the determinant of the coefficients of A; in 

Eqs. (17a-b) and Di(w) is the modification in D(ui) re- 
sulting when the Zth column is replaced by (Fg, 0, 0, 0) 
for the longitudinal motion and (0, 0, .., 0, Fq, 0, 0, 0) 
for the transverse motion respectively. Ai are complex 
quantities when 7' ^ and the formalism developed 
above allows to calculate amplitudes and phases. The 
analytical expressions for A\ are reported in Appendix 
C. 

In Fig. 5(a-b) we show the amplitudes of the displace- 
ments for the longitudinal and transverse motion as a 
function of particle position along the chain. The plots 
clearly show an exponential decay. Regarding the dis- 
placements in the longitudinal modes, there are edge ef- 
fects, which disappear if the particle which is excited is 
being in the center of the chain. In principle this cannot 
be realized in the experiments, however except for the 
first two particles the amplitudes of the displacements 
have the same damping rate and the same phase angles 
both if the excited particle is the one at the end of the 
chain and the one in the middle. This is why for all the 
calculations we have considered a force FQe~ luJt Si^/2 in- 



stead of Fq e 



- iult S liN in Eq. (17a). 



Fitting the amplitudes to an exponential curve yields 
fej. In order to find the dispersion of k ri we calculated 
the wave's phase <j) as a function of the position and fit- 
ted it to a straight line. The definition of phase velocity, 
as being the ratio between the frequency and the wave 
number, indeed, yields k r a = A(f>. In Fig. 6 the phase 
angle as a function of the distance is plotted. It is inter- 
esting to observe that k r and fcj are independent of the 
intensity of the driving force Fq, as expected in a har- 
monic model, and the results do not change if instead 
of Foe~ luJt , that is a complex force, we consider a real 
force Iq sin wt, as in the experiment; what actually plays 
a central role is just the driving frequency. Optical and 
longitudinal waves both propagate away from the excita- 
tion region; they are backward and forward, respectively. 
It should be noticed that in a ID chain with finite length 
one should expect that only standing waves would be al- 
lowed, the effect of gas damping is the suppression of the 
reflected wave from the chain's end. 

It is interesting to observe that for low densities the 
two calculational methods, that is the one in which the 
driving force is neglected and the one in which it is ex- 
plicitly taken into account, give the same results for the 
dispersion curves, when only first neighbor interactions 
are included. In Fig. 7 the results of the two approaches 
for the optical branch are compared. Note that outside 
the bond defined by the two dotted horizontal lines the 
phonon mode is strongly damped. 
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FIG. 5: Amplitude of the displacements of particles as a func- 
tion of the distance from the location of the driving force: (a) 
displacements along the ^-direction when the excited parti- 
cle is at the extremity of the chain, (b) displacements along 
y direction when the excited particle is in the middle of the 
chain. Exponential fits to the numerical data are shown by 
the dotted curves. 



Another remarkable effect, which reflects the 
anisotropy of the system, is observed with increasing 
density: the profile of displacements for the longitudinal 
mode is no longer a pure damped exponential (see Fig. 
8(a)), because reflected waves from the chain's end start 
to appear, while the amplitudes for the transversal mode 
are still exponentially decaying. We can infer that the 
effect of damping is not simply due to friction, but also 
to the external confining potential. This is confirmed 
by the calculation of the amplitude profile when 7' = 0, 
which is reported in Fig. 8(b): even in the absence of 
the friction and in the case the driving frequency is 
low enough, an exponential decay of amplitudes with 
distance is still found. 

When we include a perpendicular magnetic field, a cou- 
pling is introduced between motion along x and along 
y direction. Apart from increased mathematical com- 
plexity, the scheme developed before is still valid. It is 
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FIG. 6: Phase angle as a function of the distance from the lo- 
cation of the driving force. The calculation are done for a sys- 
tem of N — 29 particles. Optical and longitudinal waves both 
propagate away from the excitation region; they are backward 
and forward, respectively. 
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FIG. 7: Comparison between the standard calculation, in 
which the external driving force is neglected, and the "exact" 
one, in which the external driving force is explicitly consid- 
ered. The region on the left of the vertical dotted curve is 
the region in which the condition kid < n/h e is fulfilled; the 
region between the horizontal dotted curves is the region in 
which the sums are convergent. Note that the two approaches 
give the same results in the region of convergence. 



FIG. 8: (a) Profile of the amplitudes of displacements in the 
longitudinal mode as a function of the density. For high densi- 
ties the profile is no longer a simply decaying exponential; (b) 
Profile of the amplitudes of displacements in the transversal 
mode as a function of the density in the case without friction. 
Even in the absence of friction the profile is still a decaying 
exponential. 



possible to infer the real part of the dispersion relations 
from the phase angles and the imaginary part from the 
amplitude of displacements. In the absence of friction 
and driving force, the optical and acoustical branches 
are confined in different frequency bands (see thin solid 
curves in Figs. 9(a,b)), which do not cross each other and 
with a prohibited gap • The optical frequencies follow 
the cyclotron frequency and for very high field strength 
there is no significant difference between u> op t and uj c . 
The acoustical frequencies, on the other hand, decrease 
with magnetic field strength. The gap between the opti- 
cal branch and the acoustical one for large magnetic field 
approaches lo c . 

In the presence of friction and driving force, there are 
drastic changes in the dispersion relations. The frequen- 
cies are no longer confined in different bands, because 
the frequency of oscillation is the one of the external 
force, which can be varied continously. The results of 
the calculations for different intensities of the magnetic 
field are shown in Figs. 9(a-b). There are no significative 
differences in the behavior of the real part of the disper- 
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FIG. 9: (a) Real part of the dispersion relations for a weak 
perpendicular magnetic field; (b) Real part of the dispersion 
relations for a strong perpendicular magnetic field. The insets 
show the imaginary part of the dispersion relations. 



FIG. 10: Theoretical (curves) dispersion relations for the 
acoustical mode and comparison with experimental data 
(symbols) from Ref. [7] for the one chain structure: (a) real 
part of the acoustical branch; (b) imaginary part of the acous- 
tical branch. 



sion relations with respect to the case without magnetic 
field. The imaginary dispersion relations clearly show 
that the waves are overdamped in the band gap, region 
where large values for kid are found. Notice that friction 
reduces the slope of the acoustical branch in the small 
k r a region. When the curve enters the gap region it be- 
comes strongly damped as is clearly seen from the inset of 
Figs. 9(a,b). The optical mode is for all frequencies more 
strongly damped than the acoustical one. Furthermore, 
the dispersion of the optical branch is strongly modified 
by friction, i.e. it attains a negative dispersion for all 
values of the frequency. 



V. COMPARISON WITH EXPERIMENT 

In Figs. 10 and 11 the real and imaginary dispersion re- 
lations for the acoustical and optical modes for the single 
chain configuration are respectively presented, for differ- 
ent values of the parameters. 

The calculated dispersion relations are compared with 



experimental data of Ref. |jj . The experimental data are 
in good agreement with the theoretical calculations, al- 
though the system realized in laboratory is slightly differ- 
ent from the one investigated in the theory. As a matter 
of fact, in the experiment the interparticle spacing was 
not uniform: it was 15% smaller in the center than in the 
chain's end. Due to the strong damping this density gra- 
dient is not very important in the considered forced oscil- 
lation. As in the case without friction, the optical mode 
has negative dispersion, while the longitudinal one has 
positive dispersion. The dispersion depends on density, 
and therefore on interparticle spacing. For the acoustical 
mode the frequencies of vibrations increase with decreas- 
ing interparticle distance, while for optical vibration the 
frequencies decrease with increasing densities. Further- 
more, for low n e the exponential decay is stronger in 
both cases, which implies a highly damped wave. These 
findings can be easily explained because for smaller in- 
terparticle distance the interaction forces are larger, or 
in other words for low densities the interaction between 
the particles is rather weak and consequently the effect 
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FIG. 11: The same as Fig. 10 but now for the optical branch. 



of a local perturbation is less disruptive for the other 
particles. From Fig. 11(b) it is seen that the optical 
mode is mostly constrained to a central frequency band. 
Comparing the optical branches in the absence of friction 
(see Fig. 7 of Ref. |is|) with the one in the presence of 
gas damping, it is observed that with damping the wave 
propagates beyond the frequency band allowed in the ab- 



sence of damping. For k r a = ir, 



'-'opt 



is always equal to 



zero when 7' ^ independent of the experimental pa- 
rameters. This means that in the presence of damping, 
the softening of the optical mode does no longer signal a 
structural phase transition from a single chain structure 
to a double chain structure. 

It is, however, interesting to study the behavior of the 
dispersion relations when the density approaches the crit- 
ical value for which the continuous structural transition 
from the one chain to the two chain structure occurs. 
The results are shown in Fig. 12. The optical branch 
softens when approaching the critical density, while the 
acoustical branch is hardened. Notice that at the phase 
transition point: i) the real part of the optical dispersion 
becomes linear for k r a > 1.5, ii) there is a drastic change 
of slope in the optical imaginary dispersion, and iii) the 
optical mode becomes less damped. The real and imag- 
inary acoustical dispersions are less strongly influenced 
near the zig-zag transition. This can be easily explained 



s 

3 



S 

s 




FIG. 12: Dispersion relations for different values of the den- 
sity approaching the critical values for which the zig-zag tran- 
sition from 1 chain structure to 2 chain structure is observed: 
(a) optical branches; (b) acoustical branches. In the insets 
the imaginary parts of the dispersion relations are plotted. 



by the fact that the zig-zag transition, which is respon- 
sible of the splitting of the chain, acts in the y-direction. 
Therefore, signature of the zig-zag transition are more 
easily detected in the optical phonon mode. 



VI. FORCED OSCILLATIONS IN A 
MULTI-CHAIN STRUCTURE 

In Fig. 13 and Fig. 14 the real and imaginary part of 
the dispersion relations for the forced oscillations of the 
2 and 3 chain configurations are reported respectively. 
We used the approach given in the first part of Sec. IV. 
Therefore, the dispersion relations are in Fig. 13 and 14 
only given in that part of the Brillouin zone, where the 
sums in Eqs. (14a-b) are convergent. 

From Fig. 13(a) and Fig. 14(a) it is evident that 
there is a remarkable difference in the optical branches 
between the single chain and the 2 and 3 chain struc- 
tures. In the first case the optical mode has negative 
dispersion as stated before, while for the 2 and 3 chain 
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VII. CONCLUSION 



structures the optical frequencies do not exhibit a mono- 
tonic behavior. Such feature can be attributed to the 
fact that for the single chain configuration in the case 
of the transversal mode the restoring force is only due 
to the parabolic confining potential, while in the multi- 
chain configuration the restoring force depends both on 
the external confinement and on the particle repulsion. 

Fig. 13(b) and Fig. 14(b) exhibit some similarities 
with the single chain case: i) for the acoustical modes 
the damping is an increasing function of the driving fre- 
quency, ii)the optical modes are mostly constrained to 
a frequency band, and iii) the optical modes are more 
strongly damped. The approach used for the calculation 
of the dispersion relations for the multi-chain is the same 
considered in Sec. IV for an infinite number of particles, 
this is why in Figs. 13 and 14 the dispersion relations are 
not presented in whole first Brillouin zone, but only in 
that range of the frequency corresponding to convergent 
sums. 



The ground state and the normal modes of a Q1D 
multi-chain system can be studied experimentally in a 
dusty plasma, where /xm sized particles are externally 
confined by electric fields in the sheath above the lower 
electrode. The sheath conforms to the shape of the elec- 
trode, so building up an electrode with a groove-shaped 
depression in one direction allows the realization of a 
parabolic confining potential and, as a consequence, the 
formation of a chain-like crystal in that direction. 

We investigated the structural properties and the nor- 
mal modes of such a classical Q1D system of particles 
interacting through a Yukawa-type potential. The struc- 
tural transitions are of first (primarily) and second order. 
The normal modes of the system were calculated first ne- 
glecting the effects of dissipation induced by gas drag and 
then considering explicitly the presence of friction. The 
normal modes consist of longitudinal (acoustical modes) 
and transversal (optical modes). The number of acousti- 
cal branches is equal to the number of optical branches 
and is equal to number of chains in the system. In the 
presence of friction, the free oscillations of the system 
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are exponentially damped in time. The effect of a con- 
stant magnetic field on the dispersion relations was in- 
vestigated and we found that the acoustical and optical 
branches no longer cross. 

Particular attention was paid to the case of forced os- 
cillations induced by an external driving force, as was 
investigated in the experiments. We found that earlier 
approaches to calculate the phonon dispersion relations 
are no longer valid. Our theoretical results are compared 
with experimental data and a remarkably good agree- 
ment between theory and experiment was found. 

Finally, we made predictions for the single chain dis- 
persion relations in the presence of a perpendicular mag- 
netic field and for the multi-chain dispersion relations 
when the modes are excited by an external driving force. 
We found some substantial differences as well as some 



similarities in the dispersion relations between the single 
and multi-chain structures. 
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APPENDIX A 



The equations of motion for particles in the three-chain configuration in the presence of friction are in matrix form: 



/x> W \ 
•>(i) 

V'n 
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•>(2) 

V n 

x n 
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where the superscript labels the row which the particle is sitting in. The coefficients are: 

-3jK,/h 



Bi = 



B 2 



■2 2 
Jit 

2 + 6^- + 9^ 



B 4 



Be 




9«ri2 + Qn 2 r 2 12 + g 



+ 3 




+ 3 



13 



1 + 



3K7-13 



(Al) 



where T\ 2 = \] (j + 1/2) 2 + c 3 2 and r i3 = V j 2 + 4c 3 2 . 

The equations of motion for the 2 chain structure can be obtained by the 4x4 submatrixes which are included 
in the top left part of the matrix in Eq. (Al): the coefficients involved in this case are B\, B 2l B 3 and B 4l with the 
substitution h e /3 — > n e /2. 

In the presence of a constant magnetic field B = (0, 0, B) the equations of motions for the two and three chain 

• (*) • (*) 

structure are easily obtained from the equations of motion with B = 0, adding the coupling terms y' n uo' c and —x' n uj' c 
to the equations for x and y motion respectively, for particles sitting in the i th row. 
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Obviously, the case without gas drag is immediately recovered setting 7' = 0. 



APPENDIX B 

The eigenfrequencies in the 3 chain configuration are: 

w&> = VBi - B b - 7 V4 - i 7 /2, 
w£> = a/1 + B 2 - B 6 - 7 V4 - i 7 A 

= \M + B5/2 + \Jb 2 + 8B 2 /2 - 7 2 /4 - i 7 /2, 

= \fl + B* + Be/2 + y/B* + 8B 2 /2 - 7^/4 - » 7 /2, 

= i? 5 /2 - /b 2 + 8B|/2 - tV4 - H/2, 

- yJ± + B 2 + B 6 /2 -y/B% + 8B 2 /2 - 7 2 /4 - i 7 /2, 
where the coefficients B n are the same as in Appendix A. In the 2 chain configuration the eigenfrequencies are: 

w& = v/BiTB^V4- J7 /2, 

wg> = + B 4 - 7 2 /4 - *7A 

w$ = v/1 + B[ - B 3 - 7 2 /4 - i 7 /2, 

w£> = ^1 + B 2 -B 4 - 7 2/4 - i 7 /2. 

In this case the coefficients B n are obtained from the coefficients in Appendix A with the substitution n e /3 — > n e /2. 

APPENDIX C 

We present the analytical expressions for the displacements y4^' y calculated from the Kramer's rule (Eq. (42)) in 
the case of a system of N — 29 particles. The central particle at which the driving force is acting is labelled with 
N = 0. For reasons of symmetry A^f = A x ' y . 

A x ' y = (a 2 - b 2 )(-b 2 -ba + a 2 )(-b 2 + ba + a 2 )(b 4 + Ab 3 a - 4b 2 a 2 - ba 3 + a 4 )(b 4 - 4b 3 a - 4b 2 a 2 + ba 3 + a 4 )C/D; 
A x { v = -ab(-b 3 - 2b 2 a + ba 2 + a 3 (b 3 - 2b 2 a - ba 2 + a 3 )(-7b 6 + Ub 4 a 2 - 7b 2 a 4 + a 6 )C/D; 

A x ' v = ab 2 (-a 6 + a 5 b + 5a 4 b 2 - 4a 3 b 3 - 6a 2 b 4 + 3a6 5 + 6 6 )(-a 6 - a 5 b + 5a 4 b 2 + 4a 3 b 3 - 6a 2 b 4 - 3a6 5 + b 6 )C/D; 
A x,v = -a6 3 (-3fe 2 +a 2 )(a 2 -fe 2 )(-2fe 2 + a 2 )(6 4 -4fe 2 a 2 + a 4 )C/i?; 

A x ' y = -b 4 (-a 5 + a 4 b + Aa 3 b 2 - 3a 2 b 3 - 3a6 4 + b 5 )(a 5 + a 4 b - 4a 3 b 2 - 3a 2 b 3 + 3a6 4 + b 5 )C/D; 

A x ' y = -ab 5 (-a 2 -ba + b 2 )(-a 2 + ba + b 2 )(a 4 - 5a V + 5b 4 )C/D; 

A x,v = b 6 (a 2 -b 2 )(-b 3 -3b 2 a + a 3 )(b 3 - 3b 2 a + a 3 )C / D; 

A x 7 ' y = -ab 7 (-2b 2 +a 2 )(2b 4 -4b 2 a 2 +a 4 )C/D; 

Al' v = b\b 3 -2b 2 a-ba 2 + a 3 )(-b 3 -2b 2 a + ba 2 +a 3 )C/D; 

A x ' y = -ab 9 (-3b 2 +a 2 )(a 2 -b 2 )C/D; 

A x £> = b w (-b 2 + ab + a 2 )(-b 2 -ab + a 2 )C/D; 

A x [ v = ab n (2b 2 -a 2 )C/D- 

A x { y = b 12 (a 2 ~b 2 )C/D 

A x £ = -ab 13 C/D; 

Alf = -b 14 /D, 

where D = a(-3b 2 + a 2 )(5b 4 - 5b 2 a 2 + a 4 )(o 8 - 8b 6 a 2 + Ub 4 a 4 - 7b 2 a 6 + a 8 ). 
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In the case of displacements along the x-direction: a — uj' 2 + ij'w' — 2b,b= (rig/2) exp(— n/n){2 + 2n/h e + k 2 /n 2 ) and 
C = Fq, while in the case of displacements along the y-direction: a = ui' 2 — l + i"f'uj' — 2b 1 b = —(fig/2) exp(— n/h e )(l + 
n/n e ) and C — F$ . 
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